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CN \ ABSTRACT 

^ ' Low frequency detector noise in CMB experiments must be corrected to produce 

faithful maps of the temperature and polarization anisotropies. For a Planck-type ex- 
j /' i ' periment the low frequency noise corrections lead to residual stripes in the maps. Here 

I show that for a ring torus and idealised detector geometry it is possible to calculate 
^vq ' analytically the effects of destriping errors on the temperature and polarization power 

spectra. It is also possible to compute the pixel-pixel noise covariances for maps of 
C^l ' arbitrary resolution. The analytic model is compared to numerical simulations using 

a realistic detector and scanning geometries. We show that Planck polarization maps 
at 143GHz should be signal dominated on large scales. Destriping errors are the domi- 
qq , nant source of noise for the temperature and polarization power spectra at multipoles 

£ 10. A fast Monte-Carlo method for characterising noise, including destriping er- 
rors, is described that can be applied to Planck. This Monte-Carlo method can be used 
to quantify pixel-pixel noise covariances and to remove noise biases in power spectrum 
estimates. 
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1 INTRODUCTION 

This is the fourth in a series of papers (Efstathiou 2004, 2005, 2006; henceforth E04, E05, E06 respectively) outlining a 
methodology for analysing power spectra for a Planck-type cosmic microwave background (CMB) experiment*. The first 
■ of these papers discussed fast hybrid methods for estimating temperature power spectra. The hybrid method combines a 
maximum likelihood estimator at low multipoles, computed from low resolution maps, with a pseudo-C^ (PCL) estimator at 
high multipoles computed using fast spherical transforms applied to high resolution maps (see Hivon et al. 2002 and references 
therein). A near optimal hybrid estimator can be constructed if the noise properties of the input maps obey certain properties, 
e.g. if the noise is accurately white at high multipoles and if correlated noise can be accounted for in the pixel-pixel noise 
covariance matrices for the low resolution maps used in the maximum likelihood estimate. If this is the case, accurate analytic 
expressions for the PCL covariance matrices can be derived for realistic sky cuts. The hybrid estimator concept was generalised 
to polarization in E06 and analytic formulae for PCL polarization power spectra for noise-free and noisy data were derived. 
(See also the related paper by Challinor and Chon, 2005). There is a large literature on power spectrum estimation methods 
and we will not attempt to give a comprehensive bibliography here. A fairly complete set of references to previous work can 
be found in the papers of hybrid power spectrum estimators (E04, E06) . 

CMB temperature map reconstruction errors for a Planck-type were considered in E05, on the assumption that maps 
would be constructed using simple 'destriping' techniques (see e.g. Burigana et al. 1997; Delabrouille 1998; Maino et al. 1999; 
Revenu et al. 2000; Keihanen et al. 2004, 2005). Destriping algorithms are particularly suited to a Planck-type scanning 
strategy in which the sky is scanned many times on rings. In E05, it was shown that for a ring torus scanning geometry 
(and provided the instrument noise satisfies certain reasonable properties) it is possible to compute accurately the effects of 



* For an up-to-date account of the Planck satellite and its science case see 'The Scientific Programme of Planck? by the Planck Consortia 
(2005) hereafter SPP05. 
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residual striping errors on the temperature power spectrum. For realistic detector noise, the resulting error power spectrum 
is indeed accurately white at high multipoles, but at low multipoles the striping errors introduce a bias in the temperature 
power spectrum that has a distinctive 'sawtooth' type pattern with an envelope that varies approximately as l/£. For realistic 
Planck-type scanning geometries, the ring torus model provides a good (but not exact) description of the error power spectrum 
determined from simulations of the time-ordered data (TOD). E05 also argued that for a slowly precessing scanning strategy, 
and detector noise properties such as envisaged for Planck, a simple destriping map making algorithm will be close to optimal. 
The destriping errors are dominated by 'irreducable' white noise errors in the ring overlaps and so there is little to be gained 
by implementing computationally demanding least squares map making methods (e.g. Wright et al. 1996; Tegmark 1997a,b; 
Borrill et al. 2001; Natoli et al. 2001; Dore et al. 2001). This conjecture has been verified by Ashdown et al. (2006, hereafter 
A06), who compared a number of map making codes applied to simulated Planck data (specifically, simulations of four 217 
GHz detectors). The temperature power spectrum residuals from the different codes are almost identical, except for small 
differences at low multipoles. A simple destriping code (for example, the Springtide code described by A06) fitting one offset 
(or 'baseline') per ring pointing produces a map that is very close to optimal. Depending on the value of the detector noise 
'knee frequency', the accuracy of a destriping algorithm can be improved further by increasing the number of baselines and 
by applying more sophisticated destriping algorithms that include information on the correlations between baseline offsets 
(Keihanen, Kurki-Suonio and Poutanen, 2005). However, the simulations of A06 show that for the knee frequencies expected 
for Planck the improvements in the maps are minor as the number of baselines is increased above one per scanning ring. 

In this paper we extend the analysis of E05 to polarization. As we will show in Section 3, for a ring torus geometry 
it is possible to construct an analytic model for temperature and polarization destriping errors. Formulae are presented in 
Section 3 for the contribution of residual destriping errors to the temperature and polarization power spectra and also for the 
pixel-pixel noise covariances in maps of arbitrary resolution. 

Section 4 describes numerical simulations of a Planck-like focal plane geometry and slow precessing scanning strategy. 
Specifically, we use the geometry of the four polarization sensitive bolometer pairs at 143 GHz (see SPP05). However, to 
make large numbers (thousands) of simulations feasible, we have assumed 60' FWHM beams instead of the 7.1' beams of 
Planck, and that the scans consist of 1080 ring pointings instead of the 8640 ring pointings of Planck. The detector white 
noise assumed for these simulations was adjusted so that the noise properties of the maps match those expected from Planck 
at equivalent spatial resolution. These simulations therefore give a good impression of the actual noise properties of Planck 
143GHz maps at low multipoles. With these simulations it is possible to compute the pixel-pixel noise covariance matrix for 
degraded resolution maps and to isolate the component arising from destriping errors. These pixel noise covariance matrices 
can be compared with the ring torus approximation derived in Section 3. We then develop an extremely fast and accurate 
Monte-Carlo method of estimating the pixel-pixel noise covariance matrices based on the statistical properties of destriping 
baseline offsets. This method offers a reliable way of quantifying the pixel-pixel covariance matrices for temperature and 
polarization maps of arbitrary resolution and for realistic noise and scanning geometries. These pixel-pixel noise covariance 
estimates can then be used as inputs to maximum likelihood methods of power-spectrum estimation (E04, E06, Slosar and 
Seljak 2004) or map-based Gibbs sampling methods (see e.g. Larson et al. 2006). 



2 OVERVIEW OF POLARISATION MAP MAKING WITH DESTRIPING 

We follow the notation of Keihanen et al. (2005, hereafter KKP05) and write the TOD stream as 

y = Pm + n', (1) 

where m is the pixelised CMB temperature map, P is a pointing matrix, and n' is the detector noise. (Polarized maps will 
be dealt with below, but to keep the equations simple we will assume for the moment that we are solving for a temperature 
map). We assume that the noise can be decomposed into a correlated part consisting of a linear combination (described by 
the matrix F) of a set of baseline offsets a and a white noise component n, 

n' = Fa + n, (2) 
with covariance matrix 

(n'(n') T ) =FC F T + C„, C a = (aa T ), C„ = <nn T >. (3) 

As shown in KKP05, if the distributions of the offsets a and white noise n are assumed to be Gaussian, the maximimum 
likelihood map and offsets are give by minimising the x 2 > 

X 2 = (y - Fa - Pm) T C^ 1 (y Fa Pm) + a T C^a, (4) 
giving the solutions 

m = (P T C^ 1 P)- 1 P T C r ; 1 (y-Fa), (5a) 
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a = (F T C- 1 ZF + C- 1 )- 1 F T C- 1 Zy, (5b) 
where 

Z = I-P(P T C- 1 P)- 1 P T C; 1 . (5c) 

The Fisher matrix of the baseline offsets is given by the second derivative of (4) with respect to the baselines a, 

F a = (F T C^ 1 ZF + C- 1 ), (6) 

and we will assume that the posterior distribution of the offsets a is given by a Gaussian distribution of the form 

(27r|F: 1 |)- 1/2 c X p(-ia T F a a) . (7) 

The first term in (6) depends only on the scanning geometry and white noise component of the detector noise. The second 
term depends on the form of the detector noise spectrum at low frequencies. For Planck-like noise properties (detector knee 
frequencies of 30mHz), the first term is dominant and represents an irreducible contribution to the baseline offset covariance. 
In the rest of this paper we will assume that the irreducible term dominates and hence that the striping errors are independent 
of the exact form of the low frequency noise spectrum. In other words, we ignore C a and assume white noise errors in the TOD 
throughout. This should be an excellent approximation for Planck (see Figure 5 of E04). It is straighforward to generalise 
to the results of this paper include correlations C a in the baseline offsets. Methods for estimating C a from the TOD are 
discussed in KKP05. 

The generalisation of the above analysis to polarization sensitive detectors is straightforward. Define y d to be the TOD 
for detector d with orientation a d with respect to a fixed spherical polar coordinate system. In terms of the Stokes parameters 
I, Q and U in this fixed coordinate system, the detector response is 

y 



'' - \ (I + Qcos(2a d ) + Usin(2a d )) + n', (8) 



(see Figure 1 and Section 3.1 for a precise definition of the angle a). If the noise covariance C n is diagonal, C„ = o\ and the 
pointing matrix P is consists of ones and zeros (denoting when a map pixel falls onto an element of the TOD) the first term 
in equation (4) can be written as 

pd di(Zp 

where p denotes the map pixel, di C p denotes the TOD elements of detector d that lie within map pixel p and yf is the 
expected TOD signal for detector d from the best fitting maps, I p , Q p , U p , derived by minimising the x 2 of equation (4). If 
correlated errors are small, the second derivative of (9) gives the inverse covariance matrix for I p , Q p , U p , at pixel p, 



2^7% — ^ °i i 

^ cos(2a d ) -r-^ cos 2 (2a d ) ^ cos(2a d ) sin(2a ci ) 

2^ a \ ' <L — 7\ — 2^ 71 

sin(2c« d ) cos(2 Q1 ti ) sin(2a d ) sin 2 (2a d ) 



(10) 



where the sums extend over all detector TOD elements di C p. The condition number of the matrix M is often used to 
determine whether the Stokes parameters are well constrained in map pixel p (see e.g. A06). Evidently, equation (10) applies 
only if the TOD noise is strictly white and in general the noise will be correlated. Section 4 describes a fast Monte-Carlo 
method for estimating these correlated errors. 



3 ANALYSIS OF A RING TORUS 
3.1 Ring torus geometry 

In this Section, we analyse destriping errors in the idealised case of ring torus scanning geometry with two pairs of bolometers 
in the 'standard' Q and U configurations shown in Figure 1. Each detector pair is assumed to point to the same spot on the 
sky at a boresight angle of 9t to the spin axis. The sky is covered by a ring torus composed of a set of discrete N r i ng = 2ir/Aa 
rings of width Aa as the spin axis is repointed along a great circle. The detector orientations of the Q- and [/-configurations 
with respect to a fixed polarization basis on the sky can be specified by the the angles 

a q = cos _1 (— qi ■ e^,), a" = cos _1 (— ui ■ e^) — ^. (11) 

Since the detectors are assumed to point to the same spot on the sky, a q = a" for detectors in the standard configuration 
of Figure 1, though below we will write down general results for arbitrary angles. The response of each of the detectors is 
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Figure 1. Polarization orientations for two pairs of bolometers in the standard Q and U configurations as discussed in Section 3.f . The 
polarization angles a for arbitrary detector orientation are defined with respect to the unit vectors Sj, and eg defined with respect to a 
reference spherical coordinate system. 



denoted tL, tL, t* ui , t„ 2 , as shown in Figure 1. If the noise in each TOD is white, with identical variance <j\, unbiased 7, Q, 
and U maps can be constructed by minimising (9): 

i P = $3|(&+4++£i+0. (12a) 

iCp 

qPinp qp\ qP I qp i r*p\ 

Up - (SfSJ - ' 1 ' 

CP ( CP _ np\ _ qPfr'P — cp\ 

p ~ (S*S% ~ (5f) 2 ) ' 1 ' 

where 

C Lu} - ^(t I {9 , tl}l -4,, tl}2 )cos(2a^'">), (13a) 

S Lu} - ^(t l {9 , tl}l -4,, u}2 )sin(2a l { ^'" } ), (13b) 

Sf = ^(cos 2 2Qf +sin 2 2QD, (13c) 

S% = ^(sin 2 2a« +cos 2 2Qr), (13d) 

S3 = (cos 2a| sin 2af — cos 2a" sin 2a" ) . ( 13e) 

iCp 

The problem that we wish to solve in this Section is as follows: Assume that destriping is performed by fitting a single 
baseline offset for each scanning ring for each detector and assume further that the errors in the baseline offsets are uncorrelated 
and charaterised by dispersions cr 2 and <r 2 for the Q- and [/-configurations respectively. We wish to calculate the biases in the 
temperature and polarization power spectra associated with these errors. This provides a good model for the effects of residual 
striping errors on CMB power spectra. As we will show in the rest of this Section, this problem can be solved analytically for 
a ring torus geometry. 

To compare with the analytic calculations, we have run a set of 10000 simulations constructing I, Q and U maps, using 
equations (12a- 12c) for a ring torus geometry with random ring offsets. Q and U maps for one such realisation are shown in 
Figure 2. In this set of simulations, we use iVrmg = 2160 rings of width Aa = 2n/N r i ns composed of A^ing ring pixels. Maps 
were constructed using an igloo pixelisation scheme (E04) with pixels of size A8 C = 0.25°. We will refer to these maps as 
'error maps'. 

The assumptions underlying these simulations are, of course, highly idealized. As mentioned above, with these assumptions 
the error power spectra can be computed analytically. It is straightforward to generalise these simulations to realistic Planck- 
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Figure 2. Q and U maps generated from random striping errors for two polarization detector pairs in the standard configuration of 
Figure 1 and a ring torus scanning geometry with a boresight angle of 0j, = 85° . 

like scanning strategies and to include correlations between (an arbitrary number) of baseline offsets. Simulations of this sort 
are described in Section 4. As we will show, the Planck scanning strategy (see SPP05) is so close to a ring torus that the 
model described in the next Section provides quite a good description of the destriping errors expected for Planck. 

3.2 Error power spectra 

The spherical harmonic transforms of the error maps can be written as 



T 
a ( m 



^^w ik Q,i k IikYi m {9ik,<t>ik), (14a) 



ik 



ik 

a fm = ^ w ik^ik{QikRJ^ + iUikRjm), ( 14c ) 
ik 

where 

Rim = 2Yl m + -vYlm, Rgm = 2Ye m — -2Y lm . (15) 

In these equations the index k denotes the ring number, i denotes the pixel number within the ring and 0,*. is the solid angle 
of the ring pixel. Uk is the total intensity in ring k, and Qik and Uik are the Stokes parameters in ring k defined with respect 
to a fixed polarization basis. The weight factors Wik account for the averaging of the ring pixels in constructing the map and 
thus are proportional to the inverse of the 'hit count' distribution for the ring torus scanning strategy. It is useful to define 
the auxiliary functions Pe m , an< ^ G~[ m : 

Y em = AT PP (cos 9)e tm ^ = P em e im *, (16a) 
R+ m = 2V2N e ATGLe m ^G+ m e^, (16b) 
RJ m = -2^2N t A7GJ m e^ = GJ m e^, (16c) 
where 

* - (^j^f • 
» - (Iff) 1 "' 

and the functions G\ m and Gj m are 

1 ... ,s \ „m/ ^ COS t 

' • 2 

sin 



GJ m = ^ + -£(£-!)) Pr (cos 8) + (£ + m)— T7 ;Pr-i(cos 9), (18a) 



7—^— {(£- 1) cos OPP (cos 9) - (£ + m)Pel 1 (cos9)}. (18b) 
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(See, for example, Kamionkowski, Kosowsky and Stebbins, 1997, and references therein.) Note that the functions P, G + and 
G~ satisfy 

A- m = (-l) m A m , G+_ m = (-l) m G+ m , GJ_ m = -{-l) m Gi m . (19) 

To compute the power spectra of the error maps, reorient each ring to a new coordinate system (9' , 4>') in which the spin axis 
is aligned with the new z' axis. The weight function w ik is normalised so that 

Wik^lik = 4w sin 6I 6 , (20) 

ik 

hence in the new coordinate system the weight factor for each ring is given by 

Wik = i sin 9 b | sin tj>' ik \. (21) 
Define the integral I(m) to be 

I{m) = J |sin0'|e-^d^>' = ^-^(l + (-l) m ), (22) 
then the contribution of a single ring k to the T, E and B mode harmonic coefficients is 

a™ = ^sm8 b I k AaI(m)PP(cos9 b ), (23a) 



aZ k = -\sm9 b AaI(rn)[QkG+ m (cos9 b ) + iU k G; m (cos8 b )], (23b) 
aZ k = jsm9 b AaI(m)[QkGi m {cos9 b )+iUkGt m (cos9 b )], (23c) 



(II) 


= 4 (<4> + <4> + <4> + <44 


(Ql) 


= <4> + <*«>> 


(Ul) 


= <4> + <4>- 


(hQk) 


= ^«4>-<4», 


(hu k ) 


= ^«4>-<4»- 



where 

A = ^(4+4+4+4). (24a) 

Ofc = 4-4> (24b) 

Uk = 4-4> (24c) 

and the quantities tq x etc are the random numbers assigned to each detector ring in the standard configuration shown in 
Figure 1. The transformation to the new coordinate system thus leads to a simple relationship between the E and B mode 
harmonic coefficients and the quantities recorded by each detector. Evidently if the baseline offsets are uncorrelated, 

(25a) 

(25b) 
(25c) 

(25d) 
(25e) 

In the simulations described in the previous Section we assumed (t^ } = (t^ } = a\ , (t^ } = (t^ 2 ) = o\, with a q — a u = 13.6mK 
(though the actual value is unimportant for this discussion). The spherical and tensorial harmonics can be rotated to an 
arbitrary spherical polar coordinate system using the Wigner matrices D l m , m (a, (3, 7), where a, f) and 7 are the Euler angles 
specifying the rotation (see e.g. Brink and Satchler 1993; Varshalovich, Moskalev and Khersonskii 1988). Thus, each of the 
ring harmonic coefficients in equations (23a)-(23c) can be rotated back to the original coordinate system and summed to 
compute the harmonic coefficients of equations (14a)-(14c): 

at m = ^2^2a' e k ml D e m , m (a k ,l3k,^ k ) = EX/ ^" 1 ' 6 ""™ Qfc d^/ m (/3 fe )e _lm7fc , (26) 

k m' k m' 

where the d £ m , m are the reduced rotation matrices. Thus to compute the error power spectra, we set ctk = 0, /3k = tv/2 and 
evaluate the integrals over 7*.. This gives 

dJ = 2(21 + 1) \ A "^E4 m ) (Pr(cos8 b )) 2 , (27a) 

m 

Gf = 4(21 + 1) AQsin2 ^E /2 ( m ) k 2 (44(cosg 6 )) 2 + ^(G7 m (cosg 6 )) 2 ] 2 , (27b) 

© 0000 RAS, MNRAS 000, 000-000 



Polarisation Destriping Errors 7 




I 

Figure 3. Destriping power spectrum errors for the ring torus geometry analysed in the text. A boresight angle of 85° has been assumed. 
The upper panel (a) shows the results for the temperature power spectrum , the central panel (b) shows results for the polarization 
E-mode power spectrum and the lower panel (c) shows results for the B-mode. In each case, the (blue) points show the results from 10 4 
simulations with parameters as discussed in the text. The (red) solid lines show the analytic predictions of equations (27a)- (27c). The 
(green) dashed lines show the high i envelope ir/(2£ + 1) for T and 2tt/(2£ + 1) for E and B. 



G * = 4(21 + 1) Aasin2db J2 ^ [°l(G- lm {cose b )f + al(Gt m (cos8 b )) 2 ] 2 , (27c) 

Gf = -M2^ Aasin29b — 2 <4)) E /2 ( m ^ m ( coseb ) 6 ^( coseb )- ( 27d ) 

m 

Equation (27a) is identical to equation (19) of E05. Equations (27b-27d) generalise the results of E05 to polarization. 

Some features of these results (useful for understanding real experiments) are worth pointing out. If a\ — cr„, the error 
power spectra of the E and B modes are identical. Furthermore, the dominant terms in both power spectra are the G\ m 
terms which are similar to the P™ terms at high multipoles. Thus, we expect the error power spectra in the E and B modes 
to have about twice the amplitude of the temperature error power spectrum. The ring variances a 2 and a\ are determined 
largely by the white noise in the TOD and the scanning geometry (equation 6), and so differences in the detector noise levels 
will directly effect the destriping error contributions to the power spectra. Finally, if (t^ ) = (t% 2 ) , equation (27d) show that 
destriping errors on the TE (denoted X in this paper) power spectrum will be identically zero, even if the destriping errors 
are large. 

Figure 3 compares equations (27a)— (27c) with the average power spectra determined from 10 4 simulations as described 
in Section 3.1 (see Figure 1). (We do not plot Cf since the simulations assume identical ring variances for each detector and 
so the error power spectrum averages to zero. The dashed lines in Figure 3 show the envelopes C% oc l/(2£ + 1), valid in 
the limit £ 3> 1 for a boresight angle 0& = 7r/2 (E04). Evidently, the simulations and theory are in perfect agreement. The 
'sawtooth' pattern, with a l/l, envelope is a characteristic feature of map making in a total power experiment with a near ring 
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torus scanning geometry. This pattern can be seen in, for example, the analysis of Planck simulations described by A06. It is 
worth emphasing the point made in the introduction that destriping errors for a Planck-like experiment cannot be removed 
by applying an 'optimal' least-squares map making code, since it is not possible to remove the 'irreducible' striping errors 
fixed by the detector white noise and scanning. Furthermore, if the detector 1/f knee frequencies are comparable to, or less 
than, the spin frequency, simple destriping with a small number of baselines coefficients will be sufficient to get close to the 
irreducible striping errors (see e.g. E05, KKP05, A06). The analysis presented in this Section shows that for a Planck-likc 
experiment, the effects of residual striping errors in the T, E and B power spectra will be 'pinned' to the detector white noise 
level (since this largely fixes the variance of the baseline offsets) and will always dominate over the white noise level at low 
multipoles. 



3.3 Pixel noise covariances 

Optimal power spectrum estimators based on maps require the pixel-pixel noise covariance matrices. In the hybrid power 
spectrum approach discussed in E04 and E06, PCL estimators are applied at high multipoles and quadratic maximum 
likelihood (QML) estimators (e.g. Tegmark 1997c; Tegmark and de Oliveira-Costa 2001) are applied on degraded resolution 
maps to determine the power spectra at low multipoles. In the examples discussed in E04 and E06, noise was either ignored 
for the QML estimates, or assumed to be uncorrelated. The results of the previous Section show that in a realistic Planck- 
like experiment, destriping errors will dominate over white noise at low multipoles and hence must be quantified to obtain 
accurate estimates of the power spectra and their statistical distributions. The same is true for any map-based power spectrum 
estimator, including map-based Gibbs sampling methods (see e.g. Larson et al. 2006). 

For a Planck-size map, containing typically ~ 2 x 10 7 pixels it is not feasible either to compute or store a full pixel- 
pixel covariance matrix. However, since destriping errors are important only at low multipoles, all that is required for power 
spectrum estimation is an accurate model for the pixel-pixel covariance matrices in low resolution maps. In fact, for the ring 
torus geometery described in the previous Section, it is possible to calculate these covariance matrices analytically. 

The power spectra computed in the previous Section are, of course, invariant with respect to the orientation of the 
spherical coordinate system, whereas the Stokes parameters Q and U and their associated pixel covariance matrices are not. 
In the ring torus scanning geometry, the spin axis defines a fixed plane in the sky (the ecliptic plane for a Planck-type 
experiment). A common choise is to pick the z-axis of the spherical coordinate system to be perpendicular to this plane. We 
therefore rotate the maps from the coordinate system shown in Figure 2 to a new coordinate systems where the empty holes 
in the maps lie at the 'ecliptic' poles. 

We assume that the low resolution maps are constructed by first convolving the high resolution maps with a Gaussian of 
width 9 S and repixelising to lower resolution. This is equivalent to multiplying the rotated harmonic coefficients (14a)- (14c) 
by 

/, = exp(-i£(£ + l)0 s 2 ), (28) 
and synthesising new maps. Define the auxilliary functions 

= 5^P £ m'(co6 06)<&'m(*/2)/(m'), (29a) 

m' 

ft, = J2 d trn>(™s9 b )d e m , m (n/2)I(m'), (29b) 

m' 

I7 m = ^G7 m ,(cos0 6 )/ m , m (V2)/(m'). (29c) 

m' 

After some algebra, we can compute the pixel-pixel covariance matrices of the synthesised maps: 
{IJj) = ^CT 2 sin 2 6» 6 Aa^ ^ A im (cos6'0A 2 m(cos6l j )$f im $ £2 _ m / <1 /f 2 cos(m(0i - 0J-)), (30a) 

i 1 i 2 m 

(QiQi) = ^ sin 2 9 b Aa^2 E [ (Km^t-m - KA-m) (3U (i)G+ _ m (j) - - m (j)) + (30b) 

(p^mPiz-m - f^ im f7 2 _ m ) (Gi im (i)Gj 2 _ m (j) - Gt im (i)Gi 2 _ m (j))]f ei fe 2 cos(m(<^ - <fo)), 

(QiUi) = ^ 2 sin 2 6 b AaJ2J2 [ - { f tA-m - *7 im f7 2 _ m ) (G, + im «G7 2 - m (j) - GJ im (i)G+ 2 _ m {jj) + (30c) 

e ± e 2 m 

{rtmK-m - f7 im f£_ m ) (G+ m (i)G+_ m (j) - Gi im (i)Gi 2 _ m (j))]f e J t2 sin(m(& - 

where 9i, cf>i are the angular coordinates of pixel i. For simplicity, we have assumed that u 2 = u 2 = a 2 , thus (UiUj) = (QiQj). 
Figure 4 compares the analytic expressions of equations (30a) -(30d) with the results of numerical simulations. The 
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(a) simulations <II> (ta) analytic <D> 




(c) sdmulaticmB <QQ> (d) analytic 




(e) simulations <QU> (f) analytic <QU> 




50 100 150 300 50 100 150 200 



i 1 



Figure 4. Pixel-pixel covariance matrices for reduced resolution (A8 C = 10° pixels) maps for a ring torus geometry with a boresight 
angle 8i, = 85° . Since the covariance matrices are four-fold symmetric, we show just the lower quadrant of 202 X 202 elements. The figures 
to the left show the results of 10 4 simulations (the same set of simulations used to construct the power spectra shown in Figure 3) for 
the II covariance (panel (a)), QQ covariance (panel (b)) and QU cross-covariance (panel (c)). (The UU covariance looks identical to the 
QQ covariance.) The figures to the right show the analytic expressions of equations (30a)— (30d). 
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simulation results are based on the same set of simulations used to generate Figure 3, but degraded to a pixel size of 
A9 C = 10° after smoothing with a Gaussian of width 9 S = 8.5°. At this resolution a full map consists of only 404 pixels and 
because of the four-fold symmetry of the covariance matrices, we plot only the lower quadrant (pixel indices are ordered as 
i = j' + N'^{i'){i' — 1) where the index i' labels the polar angle 9, with i' = 1 corresponding to the ecliptic pole, f labels 
the azimuthal angle <j)' and N$(i') denotes the number of azimuthal angle bins at polar angle i'). The analytic results are in 
perfect agreement with the numerical simulations. The general behaviour of these pixel-pixel covariances is straightforward 
to understand. The pixel noise is highly correlated along the scanning rings, almost great circles at this resolution, modulated 
by the hit count distribution. The (QU) cross-covariances are considerably smaller than the (QQ) and (UU) covariances and 
could almost certainly be ignored for most purposes. 



4 MORE REALISTIC SIMULATIONS 

The model described in the previous Section is useful for gaining physical insight into the effects of destriping errors in 
CMB experiments. However, the focal plane and scanning geometry are highly idealised and the baseline offsets are assumed 
to be strictly uncorrelated. The model offers, at best, a first approximation to the effects of destriping errors in a realistic 
Planck-type experiment. In this Section, we discuss a simple, fast, and practical Monte-Carlo method for estimation pixel 
noise covariances for a more realistic experimental configuration. 

The focal plane geometry is modelled by the four Planck 143 GHz polarization sensitive bolometer pairs (see e.g. Figure 
3 of Delabrouille and Kaplan 2002) . This consists of two pairs of bolometers in the Q-configuration of Figure 1 and two pairs 
in the [/-configuration aligned along the scan direction. We assume a scanning pattern in which the spin axis scans along the 
ecliptic plane but with a slow precession of 5° sin(2<^>). The focal plane and scanning geometry therefore violate (mildly) the 
assumptions of the ring torus model described in the previous Section. 

The input CMB maps are Gaussian realisations of a spatially flat A-dominated CDM universe with the WMAP3 param- 
eters (fl b h 2 = 0.0223, tt m h 2 = 0.127, h = 0.73, n s = 0.95, r = 0.090, cr 8 = 0.74, in the notation of Spergel et al. 2006). For 
these simulations the B-mode anistropy is set to zero, i.e. the tensor mode amplitude is assumed to be negligible and the 
effects of gravitational lensing (see e.g. Lewis and Challinor, 2006, and references therein) are ignored. A Gaussian smoothing 
with 9 S = 0.425° was applied and maps were generated with a pixel size of A9 C — 0.5°. The input power spectra for one 
such (noise free) realisation are shown in Figure 5 and the I, Q and U maps for this realisation are shown in the left hand 
panels of Figure 7. TOD were generated for each of the eight bolometers in 1080 scanning rings consisting of 1080 ring pixels. 
The scanning rings partially overlap with a width of 30' and a separation of 20'. Only Gaussian white noise was included 
in the TOD (with an rms of 54/iK per ring pixel) so the simulations are designed to compute the irreducible white noise 
contribution to the destriping errors discussed in Section 2. Destriping was performed by determining one baseline coefficient 
for each detector scanning ring. This is done by minimising for each detector 

^ 2 = 5Z X 2a 2 n ^ ~ ak ~ gi ~ + afc ' + + A I X flfc ) ( 31 ) 

p ijkk'Gp V k / 

where n p is the hit count in map pixel p and k denotes the scanning ring. The second term enforces the condition ^2 a k = 0. 
Provided A is chosen to be large enough, the solutions for the (and their errors) are independent of A. The correction terms 
q, and q~j in (31) correct for the small orientation dependent polarization contribution to the TOD (equation 8) determined 
from the Q and U maps constructed from all of the detectors. (Maps from the destriped noisy TODs were constructed by 
a straightforward generalisation of equations (12a)- (12c) to eight detectors.) The solution for the offsets ak is therefore 
iterative: offsets are determined by first minimising (31) with qi set to zero, constructing 7, Q and U maps using these offsets 
and then recomputing offsets by using the Q and U maps to compute the polarization corrections qi. The solutions converge 
very rapidly to a stable solution. The first term in the \ 2 of equation (31) is an approximation to 

X 2 = (y-Fa) T Z T C- 1 Z(y-Fa), (32) 

derived by substituting the solution (5a) back into equation (4). 

The white noise level applied to the TOD was chosen so that the final maps have a white noise amplitude in each 30' x 30' 
map pixel about equal to that expected for the 8 polarized Planck 143GHz detectors. The simulations should thus give a good 
indication of what we might expect from Planck at low multipoles. The resolution of the simulation (set by 1080 scanning 
rings, compared to the 8640 scanning rings for Planck) has been chosen so that it is possible to generate large numbers of 
simulations of the TODs on a workstation rather than a supercomputer. It is straightforward to generalise the simulations 
described in this Section to full Planck resolution and to relax the assumptions of white noise in the TOD and one baseline 
per detector ring. 

Figure 6 shows the power spectra for the realisation of Figure 5 but now including destriping and white noise as described 
above. The maps corresponding to these power spectra are shown in the right hand panels of Figure 7. The white noise level 
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I I I 

Figure 5. The (blue) points show T, X, E and B power spectra of a single Gaussian realisation of a ACDM model with WMAP3 
parameters (see text) at a resolution of SpwHM = 1°- The tensor amplitude in this model is set to zero. The (red) lines show the 
expectation values of the power spectra for this ACDM model. The input /, Q and U maps for this realisation are shown in the left hand 
panels of Figure 7. 




I I I 

Figure 6. The (blue) points show CMB power spectrum of a noisy realisation of the ACDM model shown in Figure 5, but after 
polarization destriping as described in the text. The noise parameters, focal plane geometry and scanning strategy have been chosen to 
match the Planck 143 GHz polarization sensitive bolometers, but with the noise renormalised to compensate for the degraded resolution 
of the simulations. The (red) lines show the expectation values of the power spectra for this ACDM model as plotted in Figure 5. The 
(green) lines in panels (c) and (d) show the white noise errors on the polarisation power spectra. 
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(a) T input (d) T output 




(c) U input (f) L output 




Figure 7. The pictures to the left (a)-(c) show the noise free input /, Q and U maps used to generate the power spectra plotted in Figure 
5. The pictures to the right (d)-(f) show the output maps with realistic instrument noise and polarization destriping corresponding to the 
power spectra plotted in Figure 6. Note that the Q and U maps are dominated by white noise, though some signal from low multipolcs 
can be discerned. 

for these simulations is so low that it cannot be seen in the T and X power spectra. However, the E and B mode power 
spectra are dominated by white noise at I £ 20. The E mode power spectrum is signal dominated at £ 10 and at these low 
multipoles the dominant source of noise comes from destriping errors. Since the white noise for each detector is assumed to 
be the same, the noise contributions to the E and B power spectra are almost identical. Comparing the right hand panels of 
Figure 7 with the input maps, one can see that the Q and U maps are dominated by white noise, but some large-scale CMB 
features can be discerned in these noisy maps. 

Figure 8 shows the maps of Figure 7 smoothed with a Gaussian of width 6 S = 8.5° and repixelised onto pixels of width 
A8 C = 10°. The large scale CMB features in the output Q and U maps are now clearly visible and the dominant source 
of errors are stripes aligned (almost) along great circles in the sky. If Planck works as expected (i.e. no unexpected and 
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Figure 8. The maps from Figure 7 smoothed with a Gaussian beam of FWHM 20° and pixelised into maps with 10° X 10° pixels. The 
Planck-like polarization maps (e) and (f) are dominated by signal at low multipoles (see Figures 6c and 6d), but noise from residual 
striping errors is also visible at this resolution. 

'uncorrectable' sources of low frequency noise) the polarization maps at this frequency should be signal dominated on large 
scales. 

Figure 9 shows the error power spectra determined from 4000 simulations with parameters as described above (plotted 
as the blue circles). The (green) dashed lines in the Figure show the white noise level expected for each power spectrum. 
The (red) lines show the ring-torus model of equations (27a) - (27c) scaled to provide a good fit to the numerical results. 
(A boresight angle of 86.15°, appropriate for the polarized 143GHz Planck detectors, has been used for these predictions.) 
The ring-torus model deviates slightly from the numerical results at £ £ 5, but provides a highly accurate match to the error 
power spectra at higher multipoles. It is therefore feasible to use the error model to fit and subtract the noise biases to the 
power spectrum estimates for the E and B modes at £ ^ 5. It is also worth pointing out that the error power spectra will 
increase the variances of the power spectrum estimates approximately as 
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Figure 9. The (blue) circles show the mean error power spectra for 4000 simulations of the type shown in right hand panels of Figure 
7. The (green) dashed line shows the white noise level to which these power spectra converge at high multipoles. The (red) solid lines 
shows sum of the white noise levels and the analytic predictions of equations (27a)- (27c) for a ring-torus geometry with boresight angle 
86.15°. 



{{ACe)2) ~ (27TT) (C?MB + c|)2 ' (33) 

and will introduce coupling between coefficient of different £ (see below). 

Pixel-pixel covariance matrices averaged over 4000 degraded resolution maps, such as those plotted in Figure 4, are shown 
in Figure 10. The left hand panels show the full pixel-pixel covariance matrices, which include signal, destriping errors and 
white noise. The panels to the right show the contribution to the pixel noise arising from destriping errors alone. These can 
be compared to the equivalent figures for the ring torus geometry shown in Figure 4. Apart from differences in the amplitudes 
arising from differences between the dispersion of the ring offsets, the additional structure seen in Figure 10 is caused by the 
slow precession in the scanning strategy. 

The left hand panel of Figure 11 shows the coviariance matrix for the ring offsets [did]) for one detector determined by 
averaging over 4000 simulations. The ring offset covariance matrices are dominantly diagonal, as assumed in the ring torus 
model of Section 3, but there are some low amplitude correlations particularly with the rings separated by A0 ~ n in ecliptic 
longtitude. The panel to the right shows the covariance matrix given by the inverse of the Fisher matrix computed from the 
\ 2 of equation (31): 

Fkk = = J2 P J2 l kc P + 2A > 

F kk , = /* = y y.. lirr _ -^- + 2) 

which depends only on the scanning pattern and the TOD variance a 2 . The theoretical covariance matrix (34) provides an 
excellent match to the empirical covariance matrix derived from destriping. 




© 0000 RAS, MNRAS 000, 000-000 



Polarisation Destriping Errors 15 




(c) full cora-ariance/SO (d) striping covBriance <Q<i> 




i I 



Figure 10. The panels to the left show one quadrant of the full pixel-pixel noise covariance matrices (i.e. signal+destriping errors+white 
noise) estimated from 4000 simulations of degraded resolution destriped maps such as those shown in Figures 8d - 8f. The panels to the 
right show the contribution to these covariance matrices arising from destriping errors alone. These can be compared to the ring-torus 
model plotted in Figure 4. Note that the full covariance matrices have been divided by scaling factors to render them visible: the {//} 
covariance plotted in (a) has been divided by a factor of 5 X 10 5 ; the (QQ) and (QU) covariances in (c) and (e) have each been divided 
by a factor of 20. 
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■ ' -n ' ring i 

Figure 11. Figure 11a shows the covariance matrix of the destriped ring offsets for a single bolometer determined from 4000 simulations. 
Figure lib shows the theoretical covariance matrix computed from equation (34). 



Using l/8th of the number of rings for Planck it is feasible, as shown in this Section, to analyse noise power-spectra 
and estimate the pixel-pixel noise covariance matrices of low resolution maps by brute force simulations of the TOD. At full 
Planck resolution this will become computationally challenging, if not impossible. For example, suppose we wish to analyse 
low resolution Planck maps consisting of 10 4 pixels, one would probably require more than 10 5 simulations of the TOD at each 
frequency to accurately characterise the pixel noise covariance matrices. However, Figure 11 suggests that full simulations 
of the TOD are not required. Instead, we can decompose the noise into an 'uncorrelated' component described by the 3x3 
matrices of equation (10) and a correlated component quantifying the destriping errors given by the Fisher matrix of (34). 
We will focus on modelling the destriping errors in the remainder of this Section. 

Given the Fisher matrix (34) it is straightforward, and very fast, to generate simulations of correlated baseline coefficients 
and to compute I, Q and U maps and error power spectra. These maps can then be degraded to a lower resolution after 
smoothing by a Gaussian (28) or a more complex function (for example, if one wants to minimise contamination from the 
Galaxy). An example is shown in Figure 12. In this Figure, the left hand panels are identical to the right hand panels of 
Figure 10 show the destriping contributions to the pixel covariances derived from full simulations of the TOD. The panels to 
the right show the destriping covariances from ring simulations assuming that the ring coefficients are Gaussian distributed 
with a covariance matrix F7y. The results are essentially identical to those from the simulations of the TODs. 

Monte-Carlo simulations of the baseline offsets is therefore provide an entirely feasible approach to analysing the complex 
striping noise for Planck. The map comparisons in A06 show that for realistic '1//' detector noise it is possible to get close 
to the 'irreducible' destriping errors using a small number of baseline offsets (as mentioned in the introduction, one offset per 
ring comes close to the 'optimal' map solution) and so a manageable number of baseline offsets can be used in the Monte- 
Carlo simulations. The example shown in Figure 6 shows that the _B-mode power spectrum for the 143GHz Planck channel 
should be signal dominated at multipoles I 10 and at I 20 destriping errors will be the dominant source of noise. It is 
therefore important to model the destriping noise from Planck. The amplitude of B-mode spectrum is unknown (for a recent 
review see Efstathiou and Chongchitnan, 2006) and may well be lower than the Planck detection threshold, corresponding 
to a tensor-scalar ratio r ~ 0.1 (SPP05) An accurate analysis of destriping errors will therefore be essential for the analysis 
of primordial B-mode anistropies (particularly for the possible detection of a 'reionisation bump' at £ ;$ , which cannot be 
observed in ground based experiments). 

Finally, we mention that if the CMB fluctuations are an isotropic Gaussian random field, then we expect the harmonic 
coefficients to satistfy 

(dlrndlm*) = C(5 u i8 mm i . (35) 

The simulations described above show that the destriping errors are well approximated by a Gaussian random field and, of 
course, Gaussianity is assumed in generating simulations from the covariance matrices of the ring offsets. However, the striping 
errors, although Gaussian, are evidently anisotropic since they are highly correlated along the scan direction. Failure to model 
the destriping errors could therefore have a significant effect on tests for non- Gaussianity in the Planck data. The striping 



© 0000 RAS, MNRAS 000, 000-000 



Polarisation Destriping Errors 17 



(a) ring simulations <U> 



(b) striping covuriance <H> 
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(c) ring simulations <QQ> (d) striping covBriance <QQ> 
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Figure 12. The right hand panels are identical to the right hand panels of Figure 11 and show the pixel-pixel covariance matrices arising 
from destriping errors alone. The left hand panels show the pixel covariance matrices from 10, 000 simulations generated from Gaussian 
realisations of ring offsets with the theoretical covariance matrices of equation (34). 
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errors will also introduce correlations in the harmonic coefficients. For example, for the temperature anisotropics, E05 used 
the ring torus model to show that the striping errors lead to correlations: 

(a\ m a\? m ,) = ja 2 sin 2 9 b AaS mm , ^ P £mi P < / m2 7(mi)7(m 2 )rft lim d^ 2m /(7r/2), (36) 

mi m2 

violating the statistical isotropy of equation (35). The ring torus model of Section 3 can be used to show that destriping 
errors introduce similar correlations in the polarization harmonic coefficients. In assessing non-Gaussianity, using either map 
based statistics or statistics based on harmonic coefficients, it will be important to model the effects of destriping errors. This 
can be done at full Planck resolution by generating large numbers of simulations using ring covariance matrices, rather than 
simulations of the TODs. 



5 CONCLUSIONS 

This paper has presented an analysis of destriping errors in temperature and polarization maps for a Planck-like experiment, 
generalising the results of E05 to polarization. For a simple detector geometry and ring torus scanning pattern, we have shown 
that it is possible to compute analytically the effects of striping errors on both the power spectra and pixel-pixel covariances. 

In Section 4 we have compared the ring torus model against simulations adopting the focal plane geometry of the Planck 
polarized 143GHz detectors and a slowly precessing scan strategy. The ring torus model provides a very accurate description 
of the noise power spectra at £ £ 5. The white noise in the detector TODs in the numerical simulations was chosen to match 
the expected noise levels of the Planck detectors in 30' x 30' pixels. Thus at low multipoles, the simulations should give a 
good indication of the expected performance of Planck. The simulations show that Planck polarization maps should be signal 
dominated at large scales. The i?-mode power spectrum should be signal dominated at I 10 and destriping errors will be the 
dominant source of noise in the E- and B-mode power spectra at I ^ 20. Destriping errors will therefore need to be modelled 
accurately in the analysis of the polarization power spectra from Planck. 

As discussed in E06, since the noise is expected to be white at high multipoles, near-optimal estimates of the power 
spectra and their covariances can be obstained using a hybrid power spectrum estimator. Provided one has a model for 
the pixel covariance matrices, the low multipoles can be determined by applying a maximum likelihood estimator to low 
resolution maps. Section 4 shows that a model for the pixel covariances for maps of arbitrary resolution can be constructed 
accurately using fast simulations based on the statistical properties of the destriping baseline offsets, rather than simulations 
of the full TOD. Such fast simulations will be useful for other purposes, for example, in assessing the signficance of tests for 
non-Gaussianity in the Planck maps. 

Acknowledgements: I thank members of the Cambridge Planck Analysis Centre, especially Anthony Challinor and Mark 
Ashdown, for helpful discussions. 
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